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Abstract 

The NASA Glenn Research Center is developing advanced passive thermal management technology 
to reduce the mass and improve the reliability of space fuel cell systems for the NASA exploration 
program. An analysis of a state of the art fuel cell cooling systems was done to benchmark the portion of 
a fuel cell system’s mass that is dedicated to thermal management. Additional analysis was done to 
determine the key performance targets of the advanced passive thermal management technology that 
would substantially reduce fuel cell system mass. 


Introduction 

The purpose of this study was to analyze a single passive cooling plate that would absorb the heat 
from adjacent fuel cells within a fuel cell stack system. The results of this analysis was to provide insight 
into the requirements of a cooling plate in terms of its thermal conductivity, thickness, and specific 
gravity that would result in a passive cooling system for a fuel cell stack that would provide the essential 
uniform thermal sink while having lower mass. 


Background 

The heart of a fuel cell is an electrochemical “cell” that combines a fuel and an oxidizing agent, and 
converts the chemical energy directly into electrical power, water and waste heat. The fuel cells used by 
NASA are hydrogen-oxygen fuel cells. The fuel cells likely to be used by NASA in the future are acid- 
based Proton Exchange Membrane (PEM) hydrogen-oxygen fuel cells. An illustration of this type of cell 
is shown in figure 1 . A hydrogen molecule reacts at the anode to create a pair of protons and electrons. 
The proton ion exchange membrane conducts the protons which were generated at the anode from the 
anode to the cathode. The electrons which were also generated at the anode are conducted through the 
electrical load that is connected to the fuel cell and also reach the cathode. The hydrogen protons and the 
electrons are reacted at the cathode with an oxygen atom to produce a molecule of water. An illustration 
of a “stack” of cells connected electrically in series (shown in fig. 1) is usually employed in applications. 
The energy released as a result of this reaction is called the enthalpy of formation. When the product 
formed is liquid water (as opposed to water vapor), the enthalpy of formation is also referred to as the 
Higher Heating Value (HHV) of hydrogen. When the product formed is water vapor, the enthalpy of 
formation is referred to as the Lower Heating Value (LHV) of hydrogen. When the reaction occurs at 
298 K the energy released is called the standard enthalpy of formation. The negative sign associated with 
the Gibbs free energy indicates that this reaction is spontaneous. Not all of this change in enthalpy is 
available for the generation of electricity. The maximum amount of energy available for the generation of 
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Anode Reaction H 2 — ► 2H + + 2e _ 

Cathode Reaction % 0 2 +2H + + 2e' — ► H 2 0| iq 
Overall Reaction % 0 2 + H 2 — ► H 2 0| iq 

Enthalpy of Formation A /7° 298 = -285.83 kJ/mole 

Gibbs Free Energy AG° 298 = -237.14 kJ/mole 

Figure 1. — Proton Exchange Membrane Fuel Cell. 

electricity is the change in Gibbs free energy which is the difference between the enthalpy of formation 
and the energy losses associated with the change in entropy. This is described in mathematical terms as 

AG = AH-TAS (1) 


where 

AG change in Gibbs free energy, kJ/mole of water formed 

AH enthalpy of formation, kJ/mole mole of water formed 

T reaction temperature, K 

AS change in entropy, kJ/K/mole of water formed 

Table 1 shows tabulated values for the change in Gibbs free energy and the enthalpy of formation of 
liquid water at one atmosphere of pressure as a function of temperature. 


TABLE 1.— THERMODYNAMIC PROPERTIES 


OF LIQUID WATER (REF. 1) 


Temperature, 

K 

AG, 

kJ/mole 

AH, 

kJ/mole 

TAS, 

kJ/mole 

298.15 

-237.141 

-285.830 

-48.689 

373.21 

-225.160 

-283.454 

-58.294 


It should be noted that the AH changes very little (less than 1 percent change) from room temperature 
to the boiling point of water; also the AS is negative because the entropy of the two reactant gases is 
greater than the product water. The reduction in entropy means there must be waste heat generated in the 
reaction that is released into the ambient environment (i.e. the ambient environment’s entropy increases), 
moreover the amount of waste heat per mole increases with increasing temperature. 
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The Gibbs free energy, AG is related to the fuel cell voltage by the equation 


AG = -nFE rev (0.001 kJ/W - sec) 


( 2 ) 


where 

n number of electrons per mole of water formed (2 moles of electrons per mole of water formed) 

F Faraday’s constant (96,487 A-sec/mole of electrons) 

E rev reversible cell voltage, Volt 

Similarly, an equation can relate the enthalpy of formation, AFI, to a hypothetical, thermal neutral 
voltage. 


AH = -nFE therm (0.001 kJ/W -sec) 


(3) 


where 

E t herm thermal neutral cell voltage, Volt 

The calculation of the reversible cell voltage and the thermal neutral cell voltage based on the 
thermodynamic data contained in table 1 is shown in table 2. 


TABLE 2.— REVERSIBLE CELL VOLTAGE AND 
THERMAL NEUTRAL CELL VOLTAGE 


Temperature, 

AG, 

AH, 

F 

-■^rev? 

Etherm? 

K 

kj/mole 

kj/mole 

V 

V 

298.15 

-237.141 

-285.830 

1.229 

1.481 

373.21 

-225.160 

-283.454 

1.167 

1.469 


Both the reversible cell voltage and the thermal neutral cell voltage are commonly used. It is useful to 
compare the actual fuel cell voltage to the reversible cell voltage which represents the maximum fuel cell 
voltage obtainable. The thermal neutral voltage is hypothetical because as earlier stated, not all of the 
energy from the change in enthalpy is available to do useful work (i.e., produce a voltage). The thermal 
neutral voltage is useful for calculating the thermal efficiency of the fuel cell process and the waste heat 
generated by a fuel cell. 

Several efficiencies can be calculated. A maximum theoretical energy efficiency can be calculated 


AG _ E rev 
AH E therm 


(4) 


where 

rp maximum theoretical energy efficiency 

Using values from table 2, it can be shown that the maximum theoretical energy efficiency of a fuel 
cell is about 83 percent at 298 K and falls to about 79 percent at 373 K. An actual energy efficiency can 
be calculated 


W a I a E a 
„ _ — §_ = §_a_ 

AH I f E therm 


(5) 
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where 

T| a actual energy efficiency 

W a actual electrical energy produced, kJ/mole 

I a actual fuel cell output current, Ampere 

If theoretical fuel cell output current, Ampere 

E a actual fuel cell voltage, Volt 

It is often useful to evaluate a fuel cell’s voltage performance with what is theoretically possible. A 
voltage efficiency can be calculated 


where 

T| v voltage efficiency 



( 6 ) 


Similarly, it is useful to evaluate a fuel cell’s current performance with what is theoretically possible. 
A current efficiency can be calculated. 


r li=- i 


( 7 ) 


where 

r|j current efficiency 

Equations (4), (6), and (7) can be substituted into equation (5) to yield. 


■Ha = niil v r| t 


( 8 ) 


To analyze the thermal management, a convenient calculation of the heat generated by a fuel cell is 
needed. Since all the energy that is not converted to electrical energy is dissipated as heat, the heat 
generated by a fuel cell can be calculated as, 


Q = (AH - W a )m H , o (l000 W - sec/kj) 


( 9 ) 


where 

Q 

m H 2 0 


waste heat generated by the fuel cell, Watt 

fuel consumption, moles of water produced per second 


Substituting equations (3) and (5) into equation (9) and simplifying, 


f 


Q = 


-nFE 


Therm 


+ nF 


f \ ' 
V ! f ) 


m H 2 0 


( 10 ) 


The fuel consumption is related to the theoretical fuel cell output current as shown in equation (11). 
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( 11 ) 


If 


m. 


Substituting equations (11) into equation (10) and simplifying, 

Q-(“ If E Therm Ia E a ) 


( 12 ) 


Typically the actual fuel cell output current is approximately equal to the theoretical fuel cell output 
current. Assuming the two quantities are equal to each other, equation (12) can be further simplified to, 


Q -If(-E Ther m +E a ) 

Taking equation (13) and dividing by the active area of the fuel cell 

~ = ^-(~ E Therm + E a ) = 0 = if (“ E Therm + E a ) 


(13) 


(14) 


where 

q fuel cell waste heat density, W/cm 2 

if theoretical fuel cell output current density, mA/cm 2 

A Total active area of the fuel cell stack, m 2 

Aceii Active area of a single cell in the fuel cell stack, m 2 

Figure 2 plots the fuel cell heat generation density (the heat generated per unit of cell area) versus the 
fuel cell output current density (ref. 2). NASA’s fuel cell applications typically optimize in the lower 
current density range (<400 mA/cm 2 ), so the heat generation expected from NASA’s fuel cells is 
generally <0.3 W/cm 2 . 
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Fuel Cell Active Area Cooling Analysis 

An analytical model of a typical fuel cell cooling plate is shown in figure 3. Heat is absorbed from the 
fuel cells on either side of the cooling plate. Typically a liquid coolant is circulated within the plate and 
heat is removed convectively as the coolant passes through plate and out of the fuel cell stack. A passive 
cooling plate must conduct the heat within the plane of the plate out to one or more of the edges of the 
plate so that the heat can be transferred to a heat sink external to the fuel cell stack. Figure 4 shows the 
difference between a conventional fuel cell thermal management system and a passive fuel cell thermal 
management system. The benefits of the passive approach are reductions in mass, system complexity, and 
parasitic power as well as improvements in system reliability. 
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Pumped Loop 
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Sensors 
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Figure 4. — Fuel Cell Thermal Management Systems. 
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The key to making the passive approach workable is making the cooling plates light enough yet 
highly thermally conductive so that the heat can be effectively removed and also provide each cell in the 
fuel cell stack a thermally uniform heat sink. Analytical expressions are needed to relate the thermal 
performance of a passive cooling plate to its physical characteristics. In this analysis the cooling plate is 
modeled as a series of small planar elements stacked edge to edge. The total amount of heat entering the 
first, top-most plate (from both sides) is 


AQi - qiAAj 


(15) 


where 

AQi heat transferred through the 1 st planar element, Watt 

qi heat density absorbed by the 1 st planar element, W/cm 2 

AAi surface area of the 1 st planar element, cm 2 

The heat absorbed by the first planar element must be conducted to the next lower plate. This heat can 
be approximated as 


AQ, = 


(ktwATj ) 
ALj 


(16) 


where 

k thermal Conductivity of the planar material, W/m/K 

t thickness of the planar material, cm 

w width of the planar elements, cm 

AT i in-plane temperature differential of the 1 st planar element, K 

ALi length of the 1 st planar element, cm 


Equating AQ! from equation (15) with AQi from equation (16) yields, 


(ktwATj ) 
AL X 


= qiAAj 


Solving for AT i 


_ (qiAA 1 )AL 1 
! ktw 


(17) 


(18) 


A similar expression can be written for AT 2 , AT 3 , and AT N . Assuming ALi= AL 2 =AL 3 =AL N =AL, and 
AAi= AA 2 =AA3=AA n =AA 


at 2 = 


(qjAA + q 2 AA)AL 
ktw 


(19) 


at 3 = 


(qjAA +q 2 AA +q 3 AA)AL 
ktw 


at n = 


(qj AA + q 2 AA + q 3 AA -I 1- q N AA)AL 


ktw 


( 20 ) 

( 21 ) 
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where 

AT 2 in-plane temperature differential of the 2 nd planar element, K 

AT 3 in-plane temperature differential of the 3 rd planar element, K 

AT n in-plane temperature differential of the N th planar element, K 

q 2 heat density absorbed by the 2 nd planar element, W/cni 2 

q 3 heat density absorbed by the 3 rd planar element, W/cm 2 

q N heat density absorbed by the N th planar element, W/cm 2 

Assuming that the heat density absorbed by all the planar elements is uniform, 

qi =q2 =q 3 = -" = qN =q (22) 

If equations (19), (20), (21) and all such similar expressions for the temperature differentials of all the 
individual planar elements AT! through AT N are added together and factoring out the common factors of 
q, AA, and AL 


Z AT i 


l 


(N + (N - 1)+ (N - 2) + • • • + l)qAAAL 
ktw 


(23) 


It can be recognized that 


N 


2> T i= AT 

i 

(24) 

AA = ALw 

(25) 

AL = — 
N 

(26) 


where 

AT total in-plane temperature differential across the cooling plate, K 

L total length of the cooling plate, cm 

N total number of planar elements 

Substituting equations (24), (25), and (26) into equation (23) 

AT _ (N + (N-l)+(N-2) + --- + l)qL 2 
N 2 kt 

Lastly, the infinite series shown in equation (26) can be rewritten as 

N + (N — l)+(N — 2)-l 1-1= 

Substituting equation (28) into equation (27) yields 


(27) 


( 28 ) 
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AT = 


As N— equation (29) can be rewritten as 


(N(N-l))qL 2 

2N 2 kt 


(29) 


AT = 


qlri 

2kt 


(30) 


For fuel cell applications it is ideal to have AT=0, but practically speaking, a AT< 3 °C is acceptable. 
The reasons are two-fold, first, the fuel cell chemical process kinetics are temperature sensitive. 
Nonuniformity in the temperature of the fuel cell electrodes means some areas of the electrodes are going 
to be more active than other areas on the fuel cell electrodes. The second reason is that the water 
movement within the fuel cell is affected by local variations in the water vapor pressure (which is a 
function of the water temperature). 

The fuel cell heat loads expected from NASA applications is < 0.3 W/cm 2 . The heat loads of fuel 
cells vary depending on their operating current density, but generally speaking these heat loads are small 
when compared with microelectronics heat loads which can exceed 100 W/cm 2 . NASA’s fuel cell 
applications generally optimize at low current densities because the high voltage efficiency that results 
from operating at low current densities minimizes fuel consumption and overall spacecraft mass. Using 
the constraints on AT and q, a relationship between L, k and t is plotted in figure 5. 

Figure 5 shows that for even small fuel cells, where the heat transmission distance is 5 cm or less the 
conductivity requirements are high. While the thermal conductivity requirements are reduced for thicker 
plates, thicker plates rapidly increase the mass of the fuel cell stack. The values of the cooling plate 
thicknesses that are plotted are representative of the cooling compartments currently in fuel cell stacks. 
Figure 5 shows that copper which has a thermal conductivity of 400 W/m/K would be adequate only for 
heat transmission distances of 4 cm or less unless cooling plates >2 mm were used. Diamond, which has a 
thermal conductivity of 2300 W/m/K, is the most thermally conductive material known, and it would be 
adequate for heat transmission distances up to about 8 cm. Considering these material constraints, it is 
understandable why passive thermal systems for fuel cells have not been developed. 

Aside from providing a uniform thermal environment for the fuel cell, it is also important that the 
mass of the cooling system be minimized. A metric used to evaluate the mass properties of a fuel cell 
system’s cooling subsystem is to calculate the ratio of the waste heat handled by the cooling system to the 
mass of the cooling system components. For a typical liquid-cooled fuel cell system, the mass includes 
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the cooling cells within the cell stack as well as the cooling system components outside the cell stack 
(e.g., tubing, fittings, coolant, sensors, accumulator, etc.). This ratio should be calculated at the maximum 
heat dissipation since it is the maximum heat dissipation which sizes the cooling system. A cooling 
system mass metric of 556 W/kg was calculated for a fuel cell system tested at the NASA Glenn Research 
Center. As a goal for an advanced passive thermal management system a cooling system mass metric of 
1000 W/kg was selected. This metric can be incoiporated into equation (30) in the following manner. 

The fuel cell waste heat density, q can be thought of as the fuel cell waste heat per cell applied to the 
cooling plate divided by the planar area of the cooling plate. Using the active area of a single cell in the 
fuel cell stack as the planar area of the cooling plate, 


Q plate 
A cell 


(31) 


where 

Qpiate Fuel cell waste heat per cell absorbed by a single cooling plate, watt 

A ce u Active area of a single cell in the fuel cell stack, m 2 

Substituting equation (31) into equation (30) results in 

AT = Opiate^ ( 32 ) 

2k Acell t 

Recognizing that the area times the thickness of the plate is the plate volume, equation (32) can be 
rewritten as 

AT = Opiate^ (33) 

2kV 


where 

V total volume of the cooling plate, cm 3 

The mass of the cooling plate is related to the volume of the cooling plate if an assumption of uniform 
mass density is made. 


P = 


M 

V 


Solving equation (34) for plate volume and substituting this into equation (33) results in 


AT = 


QL 2 p 

2kM 


Recognizing the ratio of Q/M as the cooling system mass metric 

QlL 2 p 


AT = 


vMy 


2k 


(34) 


(35) 


(36) 
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Figure 6. — Cooling plate heat transmission distance, 
thermal conductivity, and plate specific gravity. 


Using the 1000 W/kg as the value of the metric, and 3 °C as the AT, a relationship between L, p, and 
the average density of the cooling plate is plotted in figure 6. 

Figure 6 shows that metals such as copper (specific gravity of 8.89 and thermal conductivity of 
400 W/m/K) and aluminum (specific gravity of 2.7 and thermal conductivity of 180 W/m/K) which are 
used in cooling applications in general, are clearly unsuitable for fuel cell cooling except for very small 
fuel cells. Table 4 lists alternative high thermally conductive materials (ref. 3). These alternative materials 
are all carbon based: diamond, diamond composites, pyrolytic graphite, and graphite composites. Highly 
Oriented Pyrolytic Graphite (FIOPG) has the best combination of high thermal conductivity and low 
specific gravity. FIOPG appears to meet the criteria if the heat transmission distance is less than 7 
centimeters. 


TABLE 4.— ALTERNATE PASSIVE COOLING PLATE MATERIALS (REF. 3) 


Material 

Thermal conductivity, 
W/m-K 

Specific gravity, 
g/cm 3 

CVD diamond 

1100-1800 

3.5 

HOPG* 

1500-1700 

2.3 

Natural graphite 

500 


Cont. CF/Cu* 

400^420 

5. 3-8. 2 

Carbon/carbon* 

400 (690) 

1.9 

Gr flake/Al* 

400-600 

2.3 

Diamond/Cu 

600-1200 

5.9 

Diamond/Al 

550-600 

3.1 


In addition to these materials, heat pipes, which use phase changes to absorb and reject heat are 
capable of extremely high effective thermal conductivities, and also are candidates for passive cooling 
plates. Fleat pipes are probably the only candidate capable of passively moving the heat for large fuel cells 
(transmission distances of >10 cm). 
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Conclusions 


This analysis of a passive approach to fuel cell thermal management leads to the following 
conclusions: 

1 . Passive cooling can be done, but in general, it requires either low heat transmission distances 
(<10 cm), very high thermal conductivity (>1000 W-m/K), or very thick cooling plates. Making thick 
cooling plates is counter to NASA’s goal of reducing the fuel cell stack mass. 

2. There are alternative materials whose thermal conductivity is high enough to move the heat 
<10 cm. These materials are diamond, diamond composites or HOPG. Heat pipes are also a potential 
alternative, and probably the only viable approach for passive heat removal in large fuel cells. 

3. The mass of a passive fuel cell thermal management system is potentially about half of that of a 
conventional fuel cell cooling system. 


References 

1 . CRC Handbook of Thermophysical and Thermochemical Data, David R. Lide, Henry V. Kehiaian, 
CRC Press, 1994. 

2. Calculated from data taken from Kenneth A. Burke, “Small portable PEM Fuel Cell Systems for 
NASA Exploration Missions”, NASA/TM — 2005-213994, A1AA-2005-5680 December 2005. 

3. Zweben, Carl “New, Ultrahigh-Thermal-Conductivity Packaging Materials,” IEEE Components, 
Packaging & Manufacturing Technology Society — Santa Clara Valley Chapter, Sunnyvale, CA, 
March 16, 2005. 


NASA/TM— 2009-215426 


12 



Appendix — Symbols 

A total active area of the cooling plate, cm 2 

A ce ii active area of a single cell in the fuel cell stack, m 2 

E a actual fuel cell voltage, Volt 

E rev reversible cell voltage, Volt 

Etherm thermal neutral cell voltage, Volt 

F Faraday’s constant (96,487 A-sec/mole of electrons) 

I a actual fuel cell output current, Ampere 

If theoretical fuel cell output current, Ampere 

L total length of the cooling plate, cm 

N total number of planar elements 

Q Waste heat generated by the fuel cell, Watts 

Qpiate Fuel cell waste heat per cell absorbed by a single cooling plate, Watts 
AQ 1 Heat transferred through the 1 st planar element, Watts 

T reaction temperature, K 

V total volume of the cooling plate, cm 3 

W a actual electrical energy produced, kJ/mole 

if theoretical fuel cell output current density, mA/cm 2 

k thermal Conductivity of the planar material, W/m/K 

^HiO fuel consumption, moles of water produced per second 

n number of electrons per mole of water formed (2 moles of electrons per mole of water 

formed) 

q fuel cell waste heat density, W/cm 2 

q | heat density absorbed by the 1 st planar element, W/cm 2 

q 2 heat density absorbed by the 2 nd planar element, W/cm 2 

q 3 heat density absorbed by the 3 rd planar element, W/cm 2 

q n heat density absorbed by the n th planar element, W/cm 2 

t thickness of the planar material, cm 

w width of the planar elements, cm 

AAi surface area of the 1 st planar element, cm 2 

AA 2 surface area of the 2 nd planar element, cm 2 

AA 3 surface area of the 3 ld planar element, cm 2 

AA n surface area of the n th planar element, cm 2 

AA surface area of the each planar element, cm 2 

AG change in Gibbs free energy, kJ/mole of water formed 

AH enthalpy of formation, kJ/mole mole of water formed 

ALi length of the 1 st planar element, cm 

AL 2 length of the 2 nd planar element, cm 

AL 3 length of the 3 rd planar element, cm 
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AL n 

length of the n th planar element, cm 


AL 

length of the each planar element, cm 


AQ, 

heat transferred through the 1 st planar element, Watts 

AS 

change in entropy, kJ/K/mole of water formed 

AT 

total in-plane temperature differential across the cooling plate, K 

AT, 

in-plane temperature differential of the 1 st 

planar element, K 

at 2 

in-plane temperature differential of the 2 nd 

1 planar element, K 

at 3 

in-plane temperature differential of the 3 rd 

planar element, K 

AT n 

in-plane temperature differential of the n th 

planar element, K 

Tla 

actual energy efficiency 


Tli 

current efficiency 


Tit 

maximum theoretical energy efficiency 


Tlv 

voltage efficiency 
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